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PREFACE 


This  study  documents  the  results  of  USAFETAC  internal  branch  tasking  (IBT)  6070506,  "Fill  in 
of  Hourly  Temperature  for  Short-Term  Missing  Data."  The  USAFETAC/DNO  analyst  was 
Major  Waller  F.  Miller. 

The  original  request  (USAFETAC  Project  60705,  "Request  for  Engineering  Weather  Data 
Support")  was  to  USAFETAC/ECE  by  HQ  USAF/LEEEU  for  calculation  of  cooling  degree 
hours  (CDH).  When  ECE  discovered  that  hourly  temperature  records  for  many  of  the  locations 
that  required  CDH  were  incomplete,  they  asked  DNO  to  provide  a  way  to  interpolate  the  missing 
data.  Problems  with  the  records  ranged  from  i.soiated  missing  reports  at  various  stations  to 
stations  witli  climatological  data  (maximum  and  minimum  values)  only.  Limited-duty  stations, 
of  cour.se,  do  not  report  at  night  or  on  weekends,  and  .synoptic  stations  report  only  every  3  hours. 

The  task  of  replacing  the  missing  or  incomplete  data  was  divided  into  two  parts:  (1)  a  short-term 
temperature  interpolator  for  less  than  6  hours  of  missing  data  and  (2)  a  long-term  climatic 
technique  to  fill  in  missing  temperatures  for  periods  greater  than  6  hours.  This  report-dealing 
only  with  the  short-term  tasking-offers  a  method  for  calculating  (interpolating)  hourly 
temperatures  when  less  than  6  consecutive  hours  of  data  is  missing. 
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SHORT-TERM  HOURLY  TEMPERATURE  INTERPOLATION 


1.  INTRODUCTION 

1 .1  The  Tasking.  The  tasking  requested  a  method  for  interpolating  hourly  temperatures  for 
the  calculation  of  Cooling  Degree  Hours  (CDH).  'Femperature  data  from  relational  database 
(DB2)  tables  was  used  to  verify  the  results.  Section  2  of  this  report  identifies  the  location.s 
cho.sen.  Section  3  describes  the  methodology  of  the  linear  and  third-order  polynomial 
interpolators,  and  Section  4  describes  verification  of  the  algorithm.  For  this  verification,  data 
gaps  of  various  sizes  were  filled,  along  with  nonstandard  hours  for  .synoptic  stations,  which 
report  only  every  3  hours.  For  gaps  up  to  6  hours,  the  Root  Mean  Square  Error  (RMSE)  remains 
less  than  2°  F,  with  a  bias  less  than  0.1°  F. 

1.2  The  Programs.  Two  Statistical  Analy.sis  Systems  (SAS)  programs  were  delivered  to 
satisfy  the  IBT  6070506  requirement.  The  first  (INTER?)  uses  a  cubic  polynomial  to  interpolate 
missing  data  points.  The  second  (TFILLER)  is  used  to  test  the  accuracy  of  the  interpolator  by 
providing  bias  and  RMSE  statistics  between  observed  and  interpolated  data.  TFILLER  al.so 
converts  hourly  observations  into  3-hourly  ob.servations  when  requested  to  do  so. 

1.3  The  Algorithm.  The  algorithm  cho.sen  for  this  project  is  designed  to  interpolate  missing 
temperature  values  for  periods  up  to  6  consecutive  hours.  One  of  the  most  common  interpolation 
algorithms  is  linear;  however,  when  working  with  nonlinear  data  such  as  the  peak  of  the 
temperature  curve,  a  large  error  can  be  introduced  by  using  points  near  the  ba.se  of  the  rapid  ri.se 
in  temperature,  as  .shown  in  Figure  I .  A  quadratic  might  be  better  near  the  maximum  or 
minimum  daily  temperature,  but  it  doesn’t  handle  the  transition  from  maximum  to  minimum 
temperature.  The  cubic  spline  or  third-order  polynomial  interpolation  has  been  .supported  by 
various  authors  (Pielke,  1984,  and  Koehler,  1977)  and  handles  continuous  values  such  as 
temperature  and  pressure  accurately. 
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TEMPERATURE  FOR  SCOTT  AFB,  IL 


□  □□ACTUAL  OOO  linear  OOO  CUBIC 

INTERPOLATION 

Figure  /.  Comparison  of  Linear  and  Cubic  Spline  Interpolation.  Five  hours  of  missing  dat;i 
were  created  between  1300  and  17(X)Z.  For  the  linear  interpolation,  one  point  on  either  side  of 
tlie  missing  data  was  used.  The  cubic  spline  used  two  points  on  either  side  of  the  data  gap. 
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2.  DATA 


2.1  Temperature  data  for  the  study  was  from  DB2  tables.  The  data  is  stored  in  tenths  of 
degree  Fahrenlieit  for  various  militaiy'  and  civilian  reporting  stations.  Temperature  data  is 
normally  available  from  1959  to  1988  (except  Scott  AFB,  for  which  only  1962  to  1985  data  is 
available). 

2.2  Weather  Stations.  USAFETAC/ECE  cho.se  the  following  stations  for  verification;  they 
are  repre.sentive  of  continental,  coastal,  desert,  subtropical,  and  tropical  regions.  Clark  AB,  RP, 
was  on  the  original  list,  but  the  required  data  was  not  in  the  standard  DB2  tables. 

STATION _  WMO#  LAT _ LON  ELEV  (Meters) 


Scott  AFB,  IL 

724338 

38®  33’  N 

89°51’ W 

138 

Pease  AFB,  NH 

726055 

43®05’N 

70®  49’  W 

31 

McChord  AFB,  WA 

742060 

47®09’N 

22°  29’  W 

98 

Luke  AFB,  AZ 

722785 

33®  32’  N 

112®23’ W 

332 

Homestead  AFB,  FL 

722026 

25®  29’  N 

90°  23’  W 

2 

Mildenhall  AB,  UK 

035773 

52®  22’  N 

0°  29’  W 

10 

Ramstein  AB,  GR 

106140 

49®  26’  N 

7®  36’  E 

238 

Honolulu,  HI 

911820 

2l®2rN 

157®  56’  W 

5 

2.3  Temperature  Units.  Temperature  was  converted  from  tenths  of  degrees  to  whole  degrees 
Fahrenheit  before  the  interpolation  began.  Specific  units  are  not  needed  for  this  algorithm. 
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3.  METHODOLOGY 


3.1  Linear  interpolation  can  be  used  to  determine  tlie  value  of  missing  temperature  for  sliort 
intervals.  Linear  interpolation,  in  its  simplest  fonn,  is  the  mean  of  the  two  data  points  on  either 
side  of  the  evenly  spaced  data  gap  that  one  is  trying  to  fill.  Linear  interpolation  produce.s  good 
results  over  short  time  spans  and  for  straight  time  series.  Most  meteorological  time  series  are  not 
linear;  the  diurnal  temperature  curve,  for  example,  is  clo.se  to  sinusoidal. 

3.2  Cubic  Polynomial  Vs.  Cubic  Spline  Interpolation.  Current  literature  makes  a  strong 
ca.se  for  using  a  cubic  spline  interpolation  for  meteorological  variables,  and  the  Statistical 
Analysis  Systems  (SAS)  currently  available  at  USAFETAC  makes  this  option  easy.  However,  a 
cubic  .spline  requires  two  data  points  evenly  spaced  on  both  sides  of  the  gap  through  which  one 
needs  to  interpolate.  In  climatological  data,  gaps  can  be  of  varj'ing  length.s,  and  the  data 
necessar>'  for  a  cubic  spline  could  be  mi.ssing.  The  cubic  polynomial  has  many  of  the  advantages 
of  the  cubic  spline,  but  the  minimum  four  data  points  do  not  need  to  be  evenly  spaced.  If 
additional  data  is  available,  that  data  can  be  used  in  the  interpolation.  For  these  rea.son.s,  a  cubic 
polynomial  was  selected  for  this  project. 

3.3  Independent  Variables.  A  third-order  polynomial  or  cubic  can  be  obtained  by  using  the 
square  and  cube  of  the  variable  as  input.  In  this  way,  SAS  linear  regression  can  be  used  to  obtain 
the  coefficients  of  the  cubic.  Our  independent  variables  for  this  study  are  the  hour,  hour-.squared, 
and  hour-culxd.  Minutes  are  converted  to  a  decimal  by  dividing  by  60  and  adding  the  decimal  to 
the  curient  hour;  for  example,  1030  becomes  10.50. 

3.4  Data  Point  (Hours)  Selection.  A  regre.s5ion  is  done  on  each  day.  The  number  of  hours 
used  in  a  regre,ssion  is  controlled  by  NUMHRS,  which  can  vary  from  5  to  15.  Five  hours  (data 
point.s)  are  the  minimum  needed  to  do  a  cubic  spline  a.ssiiming  the  middle  point  is  missing.  As 
the  number  of  points  in  a  time  .series  increa.se.s,  the  error  al.so  increa.ses.  Fifteen  is  cho.sen  as  a 
compromise  between  the  need  to  fdl  large  data  gaps  and  the  increase  in  error.  Another  check  Is 
performed  on  the  data  to  be  sure  there  are  two  valid,  ''non-missing"  ob.servations  either  side  of 
the  center  point.  This  restriction  prevents  extrapolation  from  occurring.  Each  hour  of  the  day 
becomes  the  first  value  of  the  time  serie.s.  The  start  hour  is  then  shifted  to  zero  so  that  there  is  no 
problem  cro.ssing  into  another  day  and  having  a  discontinuous  function  in  lime  (e.g.,  22,  23,  0, 
1). 

3.5  Predicted  Values.  SAS  produces  a  "predicted"  value  for  each  hour  of  the  time  .series. 
Initially  only  the  center  point  is  used  to  fill  in  the  value.  To  increa.se  efficiency,  the  u.ser  controls 
the  number  of  points  that  can  be  filled  on  each  pass,  HRINC.  The  number  of  point.s  filled  cannot 
Ire  more  than  the  number  of  hours  used.  The  prognim  takes  10  CPU  .seconds  per  hour  per  year. 
To  fill  each  hour  individually  for  a  30-year  record  takes  69  CPU  minutes.  If  you  fill  3  hours  each 
pass,  the  run  time  drops  to  23  CPU  minutc.s. 
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3.6  Evaluation,  'llie  accuracy  of  the  iiuerpoiatiou  is  evaluated  using  the  bias  and  RMSE.  The 
bias  is  defined  as  the  mean  of  tlie  error  of  estimate,  as  sliown  below: 


BIAS  = 


Wo  -Tp) 

N 


The  error  of  estimate  is  the  difference  between  the  actual  observed  temperature  (Ty)  and  the 
predicted  temperature  (7*pj  from  the  cubic  polynomial.  N  is  the  number  of  observations  used.  A 
bias  near  zero  indicates  that  there  is  an  equal  number  of  predicted  values  above  and  below  the 
actual  value.  RMSE  is  calculated  as  shown  here: 

RMSE  = 


The  RMSE  gives  a  measure  of  the  size  of  the  differences  between  ob.served  and  predicted  values: 
the  .square  of  the  RMSE  is  really  a  measure  of  the  variance  between  observed  and  predicted 
(interpolated)  temperatures. 


4.  ANALYSIS,  RESULTS,  AND  VERIFICATION 


4-1  Analysis  Summary.  For  reasons  suited  in  Section  3,  cubic  polynomials  were  selected  as 
tlie  functional  form  for  shori-tenn  temperature  intenmlation.  Cubics  can  be  fitted  to  time  series 
of  tiny  lengtli  greater  than  or  eipial  to  a  minimum  iengtii.  For  purposes  of  interpolating  across  a 
data  gap,  the  minimum  length  of  the  time  series  for  a  gap  of  length  I  is  five  poini.s:  i.e.,  one 
missing  point  plus  two  points  on  either  side  of  the  missing  point.  Accordingly,  a  minimum  time 
.series  length  of  five  points  is  u.sed  throughout  this  project.  Unknowns  to  be  determined  were  the 
optimum  length  of  the  time  series  to  be  fitted  and  the  effect  of  increasing  gap  size  on 
performance  of  the  interpolator.  These  tests  can  be  done  using  an  option  in  the  INTER? 
algorithm  that  replaces  actual  data  with  a  missing  indicator.  The  amount  of  mi.ssing  data  is  .set  by 
I IRINC  and  centered  in  the  time  series. 

4.2  Effects  of  Time  Series  Length  and  Gap  Length  on  Performance.  In  detennining 
the  most  .satisfactory  short-term  temperature  interpolator,  three  principal  efforts  were  made. 
First,  tests  were  conducted  to  determine  the  ma.\imum  length  of  time  .series  to  u.se  in  detennining 
the  cubic  polynomial  (see  4.2. 1).  Second,  with  the  ma.vimum  time  .scries  length  detennined,  tests 
were  conducted  over  a  range  of  lengths  (up  to  the  ma.ximum  detemiined  in  the  first  tc.st) 
regarding  the  effect  on  interpolation  accuracy  of  varying  the  data  gap  from  a  minimum  of  I  hour 
missing  (five  points,  including  the  gap)  to  a  ma.Kimum  of  1 1  hours  missing  (1 5  points,  including 
the  gap)  (.see  4.2.2).  On  the  basis  of  the  varying  gap  length  tests,  the  lengtli  of  t’-.e  gap  was 
limited  to  6  missing  hours.  Third,  a  special  study  was  done  of  the  6-missing-hour  case,  in  which 
performance  of  the  6-hour  gap  interpolator  was  .studied  as  a  function  of  the  length  of  the  cubic, 
from  a  minimum  of  1 0  points  to  a  maximum  of  IS,  including  missing  points  (sec  4.2.3). 

4.2.1  Tlie  first  set  of  tests  are  run  to  determine  how  long  a  lime  scries  can  be  u.sed  in  a  cubic 
polynomial  and  still  predict  temperature  accuRitcly.  In  the,sc  tests,  thi.,  number  of  hours  in  the 
time  scries  varied  from  5  to  19.  The  center  value  was  removed  and  replaced  with  an  interpolated 
value.  Tcmpcniture  data  for  Scott  AFB  for  year  1975  were  u.sed  for  this  test.  Figure  2  shows 
RMSE  as  a  function  of  the  number  of  hours  as  the  lowest  line.  When  the  number  of  hours  equals 
5,  a  true  cubic  spline  is  performed.  The  RMSE  increa.scs  with  the  number  of  hours  u.scd  and 
e.xceeds  2  for  more  than  17  hours.  The  bias  is  Ic.ss  than  a  thousandth  of  a  degree  Fahrenheit  for 
time  .scries  under  IS  hours.  Based  on  the  above  results,  we  limited  the  length  of  data  .submitted 
to  the  cubic  to  IS  hours. 
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ROOT  MEAN  SQUARE  ERROR  FOR  SCOTT  AFB,  IL 


□OO  SINGLE  OOO  maximum  OOO  SIX 

MISSING  POINTS 

Figure  2.  Root  Mean  Square  Error  as  a  Function  of  Time-Series  Length  and  Number  of 
Missing  Hourly  Temperatures  at  Scott  AFB. 

4.2.2  In  ihe  .second  sei  of  te.sis,  fhc  uppcmK>5t  line  (marked  by  dianwnds  in  Rgure  2)  i,s  il;c 
KMSE  iKsing  the  .same  data,  but  removing  all  but  four  poinLs—two  at  each  end  of  the  lime  serie.s. 
For  e.xample,  with  .i  data  points,  I  hour  i.s  .set  to  missing:  a  ma.xin’."ni  of  1 1  hours  are  .set  to 
mi.s.sing  for  15  data  points.  The  RMSE  quickly  ri.sc.s  from  1.1(5  data  poini.s).io  4.5®F  (15  data 
poiiiLs).  Similarly,  the  bias  (not  shown  in  Figure  2)  increases  with  the  size  of  the  data  gap.  For 
more  than  1 1  missing  points,  the  bias  is  more  than  a  degree. 

4.2.3  Using  results  from  the  above  test,  the  number  of  missing  observations  was  limited  lb  si.s. 
in  Figure  2  the.se  are  shown  for  NUMHRS  greater  than  1 1  by  the  short  line  marked  wi:h  circles. 
With  6  missing  hours  in  a  15-hour  linte  series,  the  RMSE  is  2.8®F.  Tire  bias  remains  a  few  tenths 
of  degrees,  with  ma.ximuni  interpolation  error  of  24'’F. 

4.3  Test  of  Interpolator  Performance  on  3-Hourly  Synoptic  Data.  To  test  the 
interpolation  technique,  TFILLER  was  coded  in  SAS.  The  purpose  of  TFILLER  is  to  convert 
hourly  data  to  3-hourIy  data  and  verify  it.  When  the  user  requests  the  3-houriy  option 
(ll  IREEHR  =  I),  the  temperturc  for  any  hour  not  divisible  by  three  is  set  to  missing.  This  new 
lempeniiiire  array  is  sent  to  the  interpolation  program,  and  a  predicted  array  is  returned.  Tlie 
predicted  ob-servaiions  from  INTERP  is  then  compared  to  the  actual  values.  Tliebias  and  RMSE 
were  calciilaicd  using  the  difference  between  actual  and  predicted  temperature. 
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4.3.1  As  mentioned  earlier,  completing  missing  hourly  observations  from  a  synoptic  station 
would  be  the  more  common  use  for  this  program.  To  test  its  accuracy,  3-hourly  data  was  createil 
from  hourly  data  for  the  eight  stations  listed  in  Section  2.  The  stations  were  chosen  to  be 
representative  of  different  climatic  regimes  with  significant  military  operations.  The  data  was 
processed  in  a  time  series  of  1 3  hours,  with  3  hours  being  predicted  during  each  regression  even 
if  a  fill  was  not  reciuired.  We  chose  13  hours  as  the  length  of  t!  time  series  because  it  would 
include  5  actual  observations  (00,  03,  06,  09,  ai:d  12Z).  If  one  of  these  observations  were 
missing,  we  still  met  the  minimum  requirements  to  do  the  cubic  polynomial:  a  gap  less  than  6 
hours  and  a  minimum  of  two  data  points  at  or  either  side  of  the  center.  Ten  is  the  minimum 
number  of  hours  required  in  the  time  series  to  have  four  3-hourly  observations. 

4.3.2  The  table  below  shows  the  results  of  this  test.  RMSE  is  between  1  and  2  for  the  entire 
data  set.  In  reviewing  monthly  RMSE,  a  few  isolated  values  less  than  1  or  greater  than  2  occur. 
Luke  AFB,  in  the  desert,  has  the  poorest  overall  RMSE.  The  bias  for  all  the  stations  still  remains 
small,  at  less  than  a  tenth  of  a  degree  Fahrenheit.  One  percent  of  the  errors  between  observed 
and  interpolated  values  were  greater  than  the  99th  percentile  and  1  percent  were  less  than  the  1st 
percentile,  as  shown.  The  maximum  error  between  the  observed  and  interpolated  for  any  hour  is 
17.3°F  at  Luke,  but  the  maximum  error  for  other  stations  is  in  the  low  teens.  With  13  hours  in 
the  time  series,  2  hours  will  be  filled  during  each  pass  for  a  synoptic  station  using  five  equally 
spaced  datii  points  for  the  regression.  A  gap  of  5  hours  may  occasionally  occur  in  the  data  when 
one  synoptic  observation  is  missing.  If  more  than  one  synoptic  observation  is  missing,  however, 
an  interplations  cannot  be  performed. 


Comparison  of  Interpolated  Temperature  using  3-Hourly  Data  to  Actual  Observations 
with  a  2-Hour  (Jap. 


Block  Percentile 


Station  Name 

Number 

Bias 

RMSE 

99 

1 

Scott  AFB,  IL 

724338 

0.002 

1.44 

3.64 

-4.07 

Luke  AFB,  AZ 

722785 

-0.006 

1.99 

4.67 

-5.58 

McChord  AFB,  WA 

742060 

-0.017 

1.31 

3.33 

-3.70 

Pease  AFB,  NH 

726055 

-0.002 

1.27 

3.-33 

-3.4 1 

MildenhallAB,  UK 

03-5773 

-0.004 

1.20 

3.01 

-3.22 

Ramstein  AB,  GR 

106140 

-0.010 

1.44 

3.70 

-4.05 

Homestead  AFB,  FL 

722026 

0.038 

1..36 

3.59 

-3.85 

Honolulu,  HI 

911820 

-0.072 

1.20 

3.07 

-3.40 

4.4  Accuracy.  The  AWS  Capabilities  Master  Plan  (AWS/CMP-89/()01)  lists  the  accuracy  of 
the  current  Air  Force  temperature  sensor,  the  l'QM-1 1,  as  1.0°F.  Its  replacement,  the  FMQ-8,  is 
also  accurate  to  l.()°F.  The  RMSE  results  from  the  tests,  then,  are  only  slightly  larger  than  the 
accuracy  of  the  equipment. 


8 


5.  CONCLUSION 


5.1  RMSE  Error.  USAFETAC/DNO  has  developed  an  effective  temperature  interpolation 
scheme  using  a  cubic  polynomial.  When  interpolating  gaps  of  less  than  6  hours,  the  RMSE  is 
below  2.8°F.  These  errors  are  randomly  distributed  between  being  over-  and  under-estimates  of 
the  temperature.  For  synoptic  data,  the  RMSE  error  remains  less  than  2°F  and  has  a  low  bias. 

5.2  Limitations.  Constraints  are  placed  in  the  software  to  limit  a  data  gap  to  6  hours  and  a 
time  series  to  1 5  hours.  These  restrictions  are  chosen  to  prevent  large  interpolation  enors.  There 
must  be  at  least  two  observations  on  each  side  of  the  middle  point  to  keep  extrapolation  errors 
from  occurring.  An  option  in  the  SAS  program,  TFILLER,  allows  the  user  to  test  the  accuracy  of 
the  interpolation  scheme  for  stations  with  hourly  observations. 

5.3  Applicability.  The  interpolator  should  be  applicable  to  any  other  continuous  atmo.spheric 
variable  such  as  pressure,  dew-point  temperature,  or  density  altitude.  The  user  should  run  a  test 
similar  to  the  one  described  in  this  report  before  proceeding;  the  software  is  designed  to  do  this 
for  any  variable  in  the  DB2  tables. 
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ACRINAB 

CDH 

CPU 

0B2 

FMQ-8 

HRINC 

INTERP 

NUMHRS 

RMSE 

SAS 

TFILLER 

TMQ-11 

Z 


SPECIALIZED  TERMS  AND  ACRINABs 

acronym,  inititalisni,  abbreviation 

cooling  degree  hour 

computer  processing  unit 

IBM  relational  database  in  use  at  USAFETAC 

AWS  temperature/dew  point  set 

The  number  of  hours  filled  during  each  step. 

SAS  .software  that  perfonns  a  cubic  polynomial  interpolation. 

The  number  of  hours  in  a  time  .series. 

Root  Mean  Square  Error 
Statistical  Analysis  Systems 

SAS  software  that  converts  an  hourly  reporting  station  to  a  3-hourly  reporting 
station  and  performs  verification. 

AWS  humidily/temperature  measuring  set 
Zulu  time 


DiSTRIBUTION 


AWS/DO.  Scott  AFB,  IL  62225-5(X)8 . 2 

AWS/XT,  Scott  AFB,  IL  62225-5008 . 2 

OL  A,  HQ  AVVS,  Buckley  ANG  Base,  Aurora,  CO  8001 1-9599 . 1 

SD/eWDA,  PO  Box  92960,  Los  Angeles,  CA  90{K)9-2960 . 1 

OL-K,  HQ  AWS,  NEXRAD  Opnl  Support  Facility,  1200  Wcsthciincr  Dr.,  Norman,  OK  73069 . 1 

OL-L,  HQ  AWS,  Kccsler  AFB,  MS  39534-5000 . . . 1 

OL-M,  HQ  AWS,  McClellan  AFB,  CA  9.56.52-5609 . 1 

Dot  1,  HQ  AWS,  Pentagon,  Washington,  DC  20330-6560 . 3 

Det2,  HQ  AWS,  Pentagon,  Washington,  DC  20330-.50.54 . 3 

Det  3,  HQ  AWS,  PO  Box  3430,  Onizuka  AFB,  CA  94088-3430 . 1 

Det9,  HQ  AWS,  PO  Box  12297,  Las  Vegas,  NV  89112-0297 . 1 

IWW/DN,  Hickam  AFB,  HI  968.53-.50(X) . 3 

1 1 WS/DON,  Elmenilorf  AFB,  AK  99.506-.50(X) . I 

20WS/DON,  APO  San  Francisco  96328-5000 . 1 

30WS/DON,  APO  San  Francisco  96301-0420 . 1 

2WW/DN,  APO  New  York  09094-5(XX) . 3 

7WS/DON,  APO  New  York  09403-.50(X) . I 

28W,S/DON.  APO  New  York  09127-50(X) . I 

3 1  W.S/DON,  APO  New  York  09 136-50(X) . 1 

3WW/DN,  OrCutt  AFB,  NE  681 13-.5(KX) . 6 

9WS/DON,  March  AFB,  CA  925i8-50(X) . . . 1 

24WS/DON,  Randolph  AFB,  TX  781.50-50(X) . I 

26WS/DON,  Barksdale  AFB,  LA  71 1 10-5002 . 1 

4WW/DN,  Peterson  AFB,  CO  80914-.5(XK) . 3 

2WS/DON.  Andrews  AFB,  MD  203.34-50(X) . 20 

.5WW/DN,  Langley  AFB,  VA  23665-.5(XX) . 3 

1  WS/DON,  MaeDill  AFB,  FL  33608-.5(XK) . I 

3WS/DON,  Shaw  AFB,  SC  29L52-5000 . 1 

.5WS/DON,  Ft  McPherson,  GA  30330-5000 . 1 

25WS/DON,  Bergstrom  AFB,  TX  78743-.5(XX) . I 

AFGWC7SDSL,  Oriutt  AFB,  NE  68 1 13-50(X) . 3 

USAFETAC,  Scott  AFB,  IL  62225-.5438 . 6 

7WW/DN.  Scott  AFB,  IL62225-.5008 . 3 

6WS/DON,  iluilburt  Field, FL 32,544-.5000 . 1 

1.5WS/DON,  McGuire  AFB,  NJ  0864I-5(X)2 . 1 

1 7 WS/DON,  Travis  AFB,  CA  94535-5986 . . . I 

3350  TCHTG/3TGU-W,  Stop  62,  Chaniite  AFB,  IL  61868-5(X)0 . 2 

3395  TCI  ITGTTrKO,  Kccsler  AFB,  MS  39534-.5(XK) . 2 

AFIT/CIR,  Wright-Patterson  AFB,  OH  4.54.33-6583 . . 1 

AFCSA/SAGW,  Washington,  DC  20.330-.50(X) . 1 

NAVOCEANCOMDET,  Federal  Building,  Asheville,  NC  28801-2723 . 1 

NAVOCEANCOMDET,  Patuxent  River  NAS,  MD  20670-5103 . 1 

NAVOCEANCOMFAC,  Stcnni.s  Space  Ctr,  MS  39529-.5(X)2  (NO  DSs,  dim  sums,  etc!!!) . I 

NAVOCEANCOMFAC,  NAS  North  Island,  San  Diego,  CA  921.35-5 130 . 1 

COMNA  VOCEANCOM,  Code  N3 1 2,  Stennis  Space  Ctr,  MS  39.529-5(XX) . I 

COMNAVOCEANCOM  (Capt  Brown,  C(k1c  N332),  Stcnni.s  Space  Ctr,  MS  39529-5(M)l . I 

N  A  VOCEANO.  Code  9220  (Tony  Ortolano),  Stennis  Space  Ctr,  MS  39.529-5(X)  I  (NO  DSs,  E3  C.! ! )  . I 

Maury  Occiinographic  Library  (NOC),  Code  XJL,  Stennis  Space  Ctr,  MS  39529-5(X)l . I 

FLENHMOCEANCEN,  Monterey,  CA  9394.3-5006 . 1 

NOARL  West,  Monterey,  CA  9.3943-5(X)6 . 1 
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